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1 Introduction 

One of the main developments in the quantum optics has been the study 
of process, for example spontaneous emission, that take place in inside, or 
adjacent to, material bodies. The need to interpret the experimental results 
has stimulated attempts to quantize the electromagnetic field (EM) in a 
material of general property [1]. 
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In an inhomogeneous nondissipative medium, quantization have been 
achieved by Glauber and Lewenstein [2]. But for dissipative medium, pres- 
ence of absorbtion has the effect of couphng the EM field to a reservoir. 
This matter is common in different quantization schemes [3|, IH O |6] since in 
contrast to classical case, losses in quantum mechanics imply a coupling to 
a reservoir whose degrees of freedom have to be added to the Hamiltonian. 
But in different methods the coupling between reservoir and EM field has 
been considered in different ways. 

In damped polarization model the EM field is coupled with matter field 
and matter field is coupled with a reservoir. This model is based on a La- 
grangian and equal-time commutation relations (ETCR) can be written be- 
tween canonical components. In this approach the Hamiltonian is obtained 
from the Lagrangian and is diagonalized by the Fano technique 

In phenomenological method the reservoir act as a noise source and these 
are conveniently represented by the Langevin force which act on the EM 
field. The EM field is written in terms of the noise operators using the Green 
function of Maxwell equations. The equivalence of these two methods has 
been shown in different papers [H El E] • 

Recently a new method for dealing with quantum dissipative systems has 
been introduced in [H [8] and has been extended to EM field in a dissipative 
media. This method is based on a Hamiltonian and the reservoir couples 
directly to EM field. In this method the time revolution of Em field, in 
contrast to damped polarization model, is obtained by solving the Heisen- 
berg equations. The extension of this method to a magnetizable media is 
straightforward [9]. 

The connection between this recent method and the other methods is the 
subject of the present paper. Throughout the paper the SI units are used. 

2 A brief review of different quantization meth- 
ods 

In this section we review the main approaches to EM field quantization in a 
dielectric medium and then compare them. 
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2.1 The damped polarization model 

This model is based on the Hopfield model of a dielectric [12]. The matter 
is represented by a harmonic polarization field and a coupling between the 
polarization field and another harmonic field is considered. Following the 
standard approach in quantum electrodynamics, quantization procedure is 
started from a Lagrangian density in real space 

^em ^mat ^r-es ^inti (l) 

where 

= - (2) 

is the electromagnetic part which can be expressed in terms of the vector 
potential A and a scalar potential U (B = VxA,E = —A — Vf7). The 
Lagrangian of the matter (dielectric) is defined by 

C-mat = — ^^X^, (3) 

where the field X is the polarization part, modeled by a harmonic oscillator 
of frequency uq. The Lagrangian 

= fM^Yl -P^Yl), (4) 

is the reservoir part, modeled by a set of harmonic oscillators used to model 
the losses. The interaction part 

/■oo 

dnt = -a(A ■ X + [/V • X) - / duv{u)X ■ Y^, (5) 

Jo 

includes the interaction between the light and the polarization field with 
coupling constant a and the interaction between the polarization field and 
the other oscillators with a frequency dependent coupling v{uj). Taking X-Y^ 
as the dissipating term is not essential but simplifies the calculations. The 
displacement field D(r,t) is given by the following combination of electric 
field and the material polarization 

D(r,t) = eoE(r,t) -aX(r,t). (6) 
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Since U does not appear in the Lagrangian, U is not a proper dynamical 
variable and the Lagrangian can be written in terms of the proper dynamical 
variables A, X and Y^^. The easiest way to do this is to go to the reciprocal 
space and write all the fields in terms of spatial Fourier transforms. For 
example the electric field can be written as 

= -At / ^^'k E(k, e^'^-^ (7) 



(2vr) 



Since E(r,t) is a real field so E (k, t) = E(— k, t). Therefore we can 
restrict integration over k to half space. The total Lagrangian is 

^ ~ J ''^^k(£em + ^mat + ^res + ^int) j (§) 

where the prime means that the integration is restricted to the half of the 
reciprocal space and the Lagrangian densities in this space are defined by 

4m = eo(E'-c'B'), (9) 

i 2 , 2 

Crnat = pX - pU^X , (10) 

Ae.= / cia;(pY^-pa;X), (H) 

JO 

tint = -o[A* ■X + A-X* + {U*X-UX*)] 

- / dujv{uj)X -Y^ + X-Y*. (12) 
Jo 

Using the Coulomb gauge k- A(k, t) = and the Euler-Lagrange equation 
for U*, we find 

t/(k,t)=z^(^^i^A^), (13) 

where 63 (k) is the unit vector in the direction of k. The matter fields X and 
Y^ can be decomposed into transverse and longitudinal parts. For example 
X can be written as 



4 



X(k,t) =x"(k,t)+X (k,t), (14) 

where k ■ X (k, t) = k x x"(k, t) = and similarly for Y^. The total 
Lagrangian can then be written as the sum of two independent parts. The 
transverse part 

where 

CL = eo{A'-c'k'A'), (16) 
= (pX^' - pa;o2X^\ (17) 

~ poo ~ I 2 

Ci,= / du{pY^' - puj'Y^ ), (18) 
Jo 

A4, = -(aA-X^*+ / (ia;t;(cu)X -Y^ + cc), (19) 
and the longitudinal part 

L\\=J'£k3, (20) 

where 



~ii ~ II ~ 112 r'^ ~ II ~ 112 

£ll = (pXll2 - p^ix" + / dupY^!^ - puWl ) 

JO 

fOO II II 

- / dcu(i;(cu)X" -Yji + cc.) (21) 
Jo 



In fl?Il) . = \JljjI + uj^, is the longitudinal frequency and a;^ = The 
link between these two parts is given by the total electric field 

E(k,t) = E^(k,t) + E"(k,t) = -A(k,t) + -X"e3(k). (22) 

Using f[T^ and the definition of the displacement field D given in ([H]), we 
recover the fact that D(r,t) is a purely transverse field as expected. 
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The vector potential in reciprocal space can be expanded in terms of the 
unit vectors ex{k) as 

A(k,t)= J2 iA(k,t)eA(k), (23) 

A=l,2 

where eA(k) • e';^(k) = 5xx' and eA(k) ■ e(k) = for A = 1, 2. 

The Lagrangian C can now be used to obtain the components of conjugate 
variables 

dC ~ 

- eo^A = ^ = eoA, (24) 
dAX 



dC 



dXl 



pXx - a Ax, (25) 



d£ ~ 

Qu^x = = pY^x - v{uS)Y^x- (26) 

'^^ujX 

Using the Lagrangian f|T5|) and the expressions for the conjugate variables 
(I24l)-(l26l). we obtain the transverse part of Hamiltonian as 

= y d^HT^im + T^mat + T^int) ^ (27) 

where 

K^^ = eo(E)^' + eoc'k'A', (28) 



is the energy density of the EM field, k is defined by A; = w A;^ + k'^ with 



k, = ^ = a/-^, and 



p Jo p 

+ rj^iiMi-'. at) +C.C., (29) 

Jo p 

is the energy density of the matter field which includes the interaction be- 
tween the magnetization and the reservoir. The frequency Uq = Uq + 
/q°° du ^^^} is the renormalized frequency of the polarization field and 
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Hi, = - j d?^A*.p + c.c], (30) 

is the interaction between the EM field and the magnetization. 

By the same method we can obtain the longitudinal part of the Hamilto- 
nian as 

~ l|2 ~ II ^ 

p Jo p 

+ r dJ-^±\\* ■ Ql) + C.C.] + l^xn (31) 
Jo p eo 

where we have separated the electric part of the Hamiltonian from the matter 
part for later convenience. 

The fields are quantized in a standard fashion by demanding ETCR be- 
tween the variables and their conjugates. For EM field components we have 

[Axik,t),El,ik',t)] = -t-S,,,Sik-k'), (32) 

^0 

and for the matter 

[X^{k,t),P*^,{k',t]=th6x,y6ik-k'), (33) 

[hxik,t),Q:,yik',t)] = ih6xx'Sik - k')Siu - uj'), (34) 

with all other equal-time commutators being zero. Indeed, it is easily shown 
that the Heisenberg equations of evolution based on the Hamiltonian fl271) and 
ETCR, (|32|) - (IM|) . are identical to Maxwell equations in a dielectric medium 
with displacement operator D(r,t), defined in 

The equations flTTI) and (jH]), together with the commutation relations in 
(!32|) — (IMI) complete the quantization procedure of transverse fields. 

In order to extract useful information about the system, we need to solve 
these equation. To facilitate the calculations, we introduce two sets of anni- 
hilation operators 

ax{k, t) = .f-^ikchik, t) + z4(k, t)), (35) 
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h{k,t) = Lf^(CooX,ik,t) + -Px{k,t)), (36) 
y 2nuJo p 

Lx{k,t) = yi£(_,^f^,(k,t) + -^Q^^{k,t)), (37) 

where (Do is defined in (!29l) . The different definitions for b and only amount 
to a change of phase and have been chosen for further simphcity. From the 
ETCR for the fields (1321) -(El we obtain 



[a,{k,t),a{ik',t)] = 6xx'S{k-k'), (38) 
[h{k,t),b{{k',t)] = 5xyS{k-k'), (39) 

[Kx{k,t),bly{k,t)] = 6xx'S{k-k')5{u - Lo'). (40) 

We emphasize that, in contrast to the previous ETCR between the con- 
jugate fields (132|) - (1M|) . which were correct only in half k space, the relations 
(138|) — (HOl) are valid in the hole reciprocal space. Inverting eq (135|) — (1371) to 
express the field operators in terms of the creation and annihilation oper- 
ators and inserting them into the matter part of the Hamiltonian ( 1271) . we 
obtain the normal ordered Hamiltonian for the transverse part of the matter 
Hamiltonian as 

4i = Tdu f d'k nwa\{uj,k)a^{u,k), (42) 
•^0 X=l,2 

Hi,, = fd'kJ2 [hCoob{{k,t)h{k,t)+ r dLohubUk,t)K,ik,t) 

A=l,2 ''° 

+ -/ du;V{u)[b{{-k,t) + b,{kMbixi-^,t) + b^,{k,t)]], 
2 Jo 

(43) 

H,nt = 4 / E miaxik) + ai(-k))(b,(k) - fel(-k)), (44) 

^ A=l,2 
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where V{uj) = ["^y^]; ^(k) = \/^^^ and integration over k has been 
extended to the whole reciprocal space. 

Now instead of solving the Heisenberg equations we use the Fano tech- 
nique to diagonalize the Hamiltonian '^L0\. After diagonalization the field 
operators are written in terms of the eigenoperators of the Hamiltonian. 

The diagonalization is down in two step. In the first step the polarization 
and reservoir part of the Hamiltonian is digitalized and the polarization is 
written in terms of its eigenoperators. Then the total Hamiltonian is again 
diagonalized with the same method. 

In the first step the diagonalized expression for H^at is obtained as (the 
calculations leading to the diagonalization are lengthy and can be found in 
[3] and we do not repeat them here) 

poo r ^ , ^ 

HLu = du d'kY: nu;Bl{u;,k)B,{u;,k), (45) 
•^0 ^ A=l,2 

where S|(k, a;) and Bx{]<.,u) are the dressed matter field creation and anni- 
hilation operators and satisfy the usual ETCR 

[Bx{uj, k, t), Bliu', k', t)] = Sxx'S{k - k')5{uj - u'). (46) 

These operators can be expressed in terms of the initial creation and and 
annihilation operators as 



roo 

B(k, u) = ao{u)b(k)+Po{up{-k)+ / dioaito, uj')b{k, uj)+p{u, uj'p{-k, to) 

Jo 

(47) 

and all the coefficient ao{u), (3o{uj), ai{u,u!') and j3i{u!,uj') can be obtained 
in terms of microscopic parameters. In [3] the relation between aQ{uj) and 
Po{u!) is obtained as 

Poiuj) = ^^^aoiu;), (48) 

U -\- UJq 

and we will use this relation in the next section. 

Using the commutators of b with B and B'' together with fj47|) . we find 

bx{k)= du[a*{Lu)Bi{k,uj)-Po{u)Bx{-k,u)]. (49) 
Jo 

The consistency of the diagonalization procedure is checked by verifying that 
the initial commutation relation between b{k) and 6^(k) are conserved. 
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Using (IHI) . (H^ and (HT]) . the total Hamiltonian can be written as 



^ = yrf3k{^kai(k)aA(k) + ^~da;;io;Si(k,u;)5A(k,Lu) 

?j /-oo 

+ -m M^?H5l(k,^)[aA(k) + ai(-k)] + //.c.} (50) 
2 Jo 

where A(A;) is defined in ( |44|) and ^^(ct;) = zao(co') + 2/?o(i^)- 

Now we should take the second step and diagonlaize the total Hamilto- 
nian, (1501) . as 

H = J d^^J hu[C^ik.uj)Cik.uj) (51) 
where eigenoperators of the system can be written as 

^(k,^) = ao{k,u)a(k) + f]o{k,uj)a\-k) 

+ / du'a{k,u,uj')B{k,u,Lo') + (3{k,u},uj')B^-k,Lj'). (52) 
Jo 

To calculate the time dependence of EM field operators we should write 
them in terms of eigenoperators of the system. 
The vector potential is given by 



A(r,t) = / d'k Y: J,^Mk,t)e^^- + H.c.]e,{k). (53) 

(27r)2 J ^^^^2 V 2eoa;k 

By inverting (152|) . writing a in terms of C and using the Hamiltonian 
( !5T|) . to calculate the time dependence of (7, we obtain d{t) as 



aA(k,t)= / du[a;{k,u)Cx{k,uj)e-''^' - Po{k,uj)Ci{k,u)e''^% (54) 
Using fl53|) and fl54|) we obtain 

A{r,t) = ( — — )2 / d'^k / du ^ 



T^W J Jo ujI + uj^{1 + x{^)) 
X C'A(k,^)e^(")e,(k) + i/.c.]. (55) 

A=l,2 
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where x{^) is obtained as 

£ >0 2 J_oo UJ — UJ — IS Z J-oo UJ ~ UJ I 

(56) 

and /(a.) = v^^(^). 

The transverse part of EM field can be obtained by the same method 
using the Hamiltonian (13T|) and the definition of E"^ in (122|) as 



E\r,t) = i^^] / duj d^k^ ^ 

X (C'3(k,u;)e*(-*-'^"-)-i/.c.)e3(k). (57) 

The polarization can be obtained by writing the matter field in terms of 
C and (7^ as 

P(r,t)= / M[eoX(^)E(r,^) + P^(r,aj)]e-*"* + /i.c.} (58) 

JO 

where 

Pnx{v,uj) = J d''\iphe^ImxC,x{K^)e'''-'- (59) 

Also using fl59|) we can show that the obtained polarizability satisfies the 
Kramers-Kronig relations. 

Although the damped polarization model is based on a microscopic model 
but as can be seen from f l55|) and ( 157|) . the final results only depend on the 
macroscopic parameter x(^)- Then this model can be used for any medium 
with known susceptibility. 

The second term of (!58|) has no classical equivalent and represents a 
Langevin fluctuation term (noise operator) which is a characteristic of a 
dissipative medium. It is easy to show that this noise operator satisfy the 
fluctuation-dissipation theorem [13]. 

2.2 Phenomenological model 

This model is based on Maxwell equations, Kubo's formula and dissipation- 
fluctuation theorem. In accordance to dissipation- fluctuation theorem, since 
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the medium is dissipative, we should add to Maxwell equations a noise field. 
This field is considered as the source of electromagnetic field. The electro- 
magnetic field can be written in terms of this noise operators by using the 
Green function of classical Maxwell equations. 

In this method the Em field operators are separated into the positive and 
negative frequency parts in the usual way. For example 

E(r,t) = ^/ dcj[E (r,cj)e-*"* + E (r,cu)e*"*], (60) 
V27r "'0 

where the positive and negative frequency parts involve only the annihilation 
and creation operators respectively. Similar decompositions hold for other 
field operators. The field operators satisfy Maxwell equations which in the 
frequency domain are of the forms 

V X E^(r, u) = za;B^(r, u), (61) 

V X B^(r, uj) = -iufiot>'^{r, u) (62) 
where the monochoromattic electric and displacement field are related through 

D^(r, u) = eo£(a;)E^(r, u) + Pl{r, uj), (63) 

the parameters eo and //q are permutabihty and permeability of free space 
respectively and e{u) = 1 + xi^j). 

In this equation, P^(r,a;) stands for noise current density operator as- 
sociated with absorption nature of the medium. The fiuctuation-dissipation 
theorem and Kubo's formula require that the noise operator satisfy the fol- 
lowing commutation relation [13] 

[P^,{r,u;),Pj,.{r',Lo')] = 2eohImx{io)6,J{r-r')6{u-uj'), 
[P+(r,^),P+^.(r,^)] = [P^,(r,a;),P^^.(r,a;)] =0. 

(64) 

In a gauge in which scaler potential vanishes, we have 

{r,uj) = iuA^ {r,u), (65) 
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B (r,a;) = V X A (r,tu). (66) 

Combining equations (!62l) . (163!) . (!65l) and (!66|) one can easily show that 
the positive frequency part of the vector potential operator satisfies 

2 

V X [V X A+(r,a;)] -eH^^A+(r,a;) = ^P;^(r,a;). (67) 

The differential equation ( 1671) can be converted to an algebraic equation 
using a Fourier transformation 

2 

- k X [k X A^(k, u)] - e{u)'^A^ik, u) = ^P;^(k, a;) (68) 

c eoC 

From ( l68l) . the positive frequency component of the vector potential in 
reciprocal space, can be written as 

A\K^)^^lLzfl^].pUk..), (69) 

where / is the unit Cartesian tensor and kk is the normal Cartesian dyadic. 
The expression in the bracket is the Green function of Maxwell equation. We 
now decompose the Longitudinal and transverse parts as 

i+-L/, s zu; — e3(k)e3(k), - + . ,„ , 

A (k,a;) = -[ , 2 2 2 r ^ ]-Piv(k,a;), (70) 

and 

A (k, u) = ■ Pjv(k, uj). 71 

Now let us define a new set of boson operators as 

MKu) = ^^ip±=. (72) 

Using the inverse Fourier transform the operators A (r, t) and A"(r,t) 
can be obtained from the equations (ITOll and ( |7T1) as 
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A (r,t) 



. h 



A=l,2 



X 



:5A(k,a;)e 



k-r) 



+ ff.c.)eA(k), 



(73) 



A"(r,t) = i^]^ Tduld^k^^ 



X (S3(k,a;)e*(-*-'^-'-) + /i.c.)e3(k), (74) 
and the electric field can be obtained from 

E(r.O^-%^. (75) 

The relations fl73|) - fl7il) are equivalent to the results of the damped polar- 
ization method (l55l) and (1571) . Therefore these two approaches are equivalent. 



2.3 Minimal coupling method 

In the minimal coupling method quantum electrodynamics in a linear polar- 
izable medium can be accomplished by modeling the medium with a quantum 
field, namely a matter field interacting with electromagnetic field. This quan- 
tum field describes the polarizability character of the medium and interacts 
with the displacement field D through a minimal coupling term. [9J (This 
method is named minimal coupling method since the matter field is coupled 
to EM field in a minimal coupling way). The Heisenberg equations for the 
electromagnetic field and matter quantum field lead to both Maxwell and 
constitutive equations. In this method we use the Coulomb gauge V ■ A = 
and the conjugate canonical momentum density of electromagnetic field is 
the displacement vector operator D(r,t) which satisfies the commutation 
relation 

[A,(r,t),Dv(r',t)] =*Mi,,(r-r'), (76) 

where Sj^^^i is the transverse delta function [H]. The total Hamiltonian is 
written as 
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^ _ / ,3 . [D(r,t)-P(r,t)p (Vx A(r,t))^ 

+ J2 [ d^^ r duhujI3i{k,u)Bx{u,k), (77) 

A=l,2,3 ° 

where Bx(k,uj) are matter field operators which satisfy the boson type com- 
mutation relations 

[Bx{'i^,Lj),Bi,{k',u;)] = Sx,x'6{uj-u;')S{k-k'). (78) 

In the Hamiltonian (I77|l . P(r,t) is the Polarization density operator of 
the medium defined by 



P(r,t)= -—J du [F(^)5A(k,a;,t)e^'^-'' + /f.c.]e,(k), (79) 

^ (27r)Wo ^^,^2,3 



where the the function F{u) is the coupling function between the electromag- 
netic field and the matter field. Electric field in terms of the displacement 
vector D and polarization vector P is defined as 

eoE(r,t) = D(r,t)-P(r,t). (80) 

Applying the Heisenberg equation to the vector potential A and using the 
commutation relations defined in ( 1761) we find that 

E^M) = %^ (81) 

and longitudinal part of electric field is obtained as 

E"(r,t)=p"(r,t). (82) 

Similarly, applying the Heisenberg equation to the displacement vector 
D we obtain 

%^ = VxH(M). (83) 
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Equation (ISUi) . (ISTI) and (155]) show that the supposed Hamiltonian (1771) gives 
the correct dynamical equation (Maxwell equation) and interpreting P as 
the polarization vector is acceptable. 

By using the commutation relation (178|) . the Heisenberg equation for 
Bx(k,(jj) can be obtained as follows 



l3x(k,Lj,t) = -iujI3x(k,uj,t) + — / dVF*(cj,r')e-'''-'"'E(r',t) •eA(k). 

(84) 

This equation have the following formal solution 

+ / TdtV^^*-*') / dr'F*{uj)e~'''-''E{r',t') ■ ex{k). 



(85) 

Using this equation and (1791) . the polarization operator can be obtained as 

P(r,t)= f'dt'x{t~t')E{r,t') + PN{r,t), (86) 
Jo 



where 

2 /■°° 



V roo 

X{t) = — duj\F{uj)\'sm{u;t), (87) 



hso Jo 

and Piv(r, t) is the noise operator of the medium and is obtained as 

r d'^\i _ ^ 

P^(r, t) = / / rfa; ^ F{uj)Bx{k, u;)e-''^'+'''-'ex{k) + //.c. 

(27!-) 2 Jo t^-, 



A=l 



From (1871) . the imaginary and real part of the Fourier transform of suscepti- 
bility can be written as 

Im[xiu;)] = ^\F{co)\\ (89) 

As can be seen from ( 189|) and (l90l) . Imx{uj) and Rex{^) satisfy the Kramers- 
Kronig relations, which verifies the interpreting of x(^) as the susceptibility. 
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In fl88l) . we can separate the positive and negative components of the 
noise operator PAr(r, t) as defined in the phenomenological method and then 
using fl89|) and flTHl) we find 



[Pmir, uj), PNj{j\ u;')] = 2eohImx{uj)5ij6{u - u')5{r - r') 



(91) 



which is the pustulated relation in the phenomenological method flMl) . 

In a physical situation where the electric susceptibility is known, the 
coupling function can be obtained from fl89l) and the Hamiltonian fl77|) can 
be constructed in terms of the electric susceptibility. Then the equation of 
motion can be obtained from flHU]) . flSTl) and fl55]) as 



In this equation if we change the lower limit of integral from to — oo and 
take the time- Fourier transform of this equation we can obtain the transverse 
part of the postulated equation of the phenomenological method fl68|) . Also 
the time- Fourier transform of the longitudinal part of electric filed, is 
equivalent to the longitudinal part of the postulated equation of the phe- 
nomenological method. Then as can be seen the phenomenological method 
can be obtained from the minimal coupling method. 

Equation (!92l) can be solved by going to the reciprocal space. The vector 
potential in reciprocal space can be written as 



where the Fourier component A can be written in terms of creation and 
annihilation operators 




at 



(92) 




(93) 



A(k,t)= Y: \ [a,(k,t)e,(k) + ai(-k,t)e,(-k)]. 




(94) 



In terms of the Fourier components, equation (!94l) is written as 
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A + ujlk+-^J^dt'xit-t')Aik,t') 
1 r°° 

= dLj[ujF{u)Bx{k,uj)e-''''ex{k) + H.c], 

(95) 



where = c|k|. 

The equation (!95l) can be solved by using the Laplace transformation, the 
details can be found in reference [9], and the final result is 



A{T,t) = fd'kJ2J—-^^ [z(a;k,t)e*^-^a,(k,0) + i/.c.]e,(k) 



x^,2 V 2(27r)3eoa;k 
1 ^ r d^k 



/fi K Z"*^ 
—-^ / du[au;,u^,t)B,{k,u,0)e''^-^ + H.c.]e,ik), 

(96) 



where 



.K,f) = L-'{4^^|^^^}, (97) 

e(^, u;k, t) = F{uj)L-'{ ' ^ }, (98) 

and L^^{.} denotes the inverse Laplace transform. The transverse electric 
field is obtained as 

E^M) = -%^, (99) 



and 



B(r,t) = V X A(r,t). (100) 
From equations ( l82l) and ( l85l) the longitudinal part of electric field is 
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E"(r,t) 



P"(r,t) 




1 

So Jo 



£o 

duo f d^k[Q{u,t)F{u)B3{k,uj,0)e'^-'' + H.c.]es{k) 



(101) 



where 



Q{io,t) = L-^{ 



1 



}■ 



(102) 



In section 4 we will show that these results, fl96|) — f fT02|) . are equivalent to 
the results obtained in the previous methods. 

3 Derivation of the minimal coupling Hamil- 
tonian from the damped polarization model 

In this section we want to obtain the Hamiltonian of minimal couplig method 
from the Lagrangian of damped polarization model. In fl32l) the conjugate 

of A is the transverse electric field — cqE but in minimal coupling scheme, 
as given in ( 1761) . the conjugate of A is the displacement vector D. We can 
change the conjugate of A from E to D by adding to the Lagrangian ([1]) 
the term a^^^^ which is a canonical transformation and do not effect the 
equations of motion. This canonical transformation leads to a A ■ X type of 
coupling which gives the displacement field — D as the conjugate of A. In 
this case using the Fourier transform the transverse part of the canonically 
transformed Lagrangian (C^) can be written as 



The other transverse or longitudinal parts of the Lagrangian do not 
change. Now C is used to obtain the components of the conjugate variables 
of the fields 




(103) 



where 




(104) 
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dC ~ ~ 
- Da = ^ = eo^A + aXi, (105) 



Pi = ^ = pXi, (106) 
dC ~ 

QL = ^ = PYJ, - v{co)Y,\ (107) 

where fll04p is used to obtaining (11051) — (11071) . Using the Lagrangian (11031) 
and the expression for the conjugate variables in (1105^ — 0107^ we obtain the 
Hamiltonian for the transverse field as 

H^ = I'd'kinL + niJ, (108) 

where 

nL = i )' + , (109) 

Co Po 

is the electromagnetic energy density. The interaction between electromag- 
netic field and the polarization field is embodied in electromagnetic energy 
density. Using ([6]) we can rewrite (I109P as 

ni-^ = eoE\k) + eoc'B\k). (110) 

This Hamiltonian is like the Hamiltonian of the free EM field. But since 
in this Hamiltonian electric field E(k) is not the conjugate component of 
vector potential A(k) then we can not use it to separate the eagenmode of 
the Hamiltonian and quantizing the EM field. 

Again the fields are quantized by demanding the ETCR between compo- 
nents and its conjugates. Then we should change the equation ( !24l) as 

[Axik,t),b*y{k',t)]=ihSxx'Sik-k'). (Ill) 

The matter and reservoir parts of the Hamiltonian do not change and we 
can use the result of section 2-1 to write 

HLu = Tdu [ d'k Y: hujBiik,u)B,ik,u), (112) 
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where ^^(k, a;) and B\{k,Ljj) are defined in (H7|) . 

The electromagnetic part of the Hamiltonian has been written in recip- 
rocal space and can be inverted to real space by using the inverse Fourier 
transform as 

J 2eo 2/io 

where X (r, t) can be written in terms of its Fourier component X"'"(k, t) as 

X^(r, t) = / d'k X^(k, t) e'^ \ (114) 

Now using the definition of h\{k, t) and &^(k, t) in we obtain X (r, f) 

as 



±\v,t) = -^J^ j rf^k 5: [&,(k,t)e^'^- + /..c.]e,(k). (115) 
(27r)2 V 2pu;o Att2 



Using (H9!) and (11 151) we have 



X^(r,t) = --^ Trfu; / d^k [-^5,(k,.;,t)e^'^-^ + if.c.]e,(k). 
(27r)2 7o J ^tt,2 « 

(116) 



where 



i^H = -\/^KM-/^oH)- (117) 

Now using ( 1116P and ( 1113^ the total transverse Hamiltonian ( llOSp can be 
written as 



rr-L rr-L i rr-L 

em "I" mat 

,D-P 



(V X A)^ 



250 2/io 
+ r duj I d^k ^ ;icj5i(k,a;)5A(k,a;,t), (118) 

A=l,2 
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where 



P fr. 



t) = / du / d'k J2 [F{u)Bx{k,u,t)e'''-' + H.c.]ex{k). (119) 

(27r)2 Jo J 



Using the same method the Longitudinal part of the matter field can be 
written as 

1 , /■ rF(uj) 

(2^ 

The total longitudinal part of the Hamiltonian ( !3TI) can be written as 



X"(r,t) = -i-^ Hduj I (i3k[:^53(k,u;,t)e^'^-' + /i.c.] e3(k). (120) 

f27rl2 Jo J a 



em. ~ ^■'■mat 



(27r)i J 2eo 
+ £ du J d^khu Bl(k,uj)B3{k,io), (121) 

where -F(co') is defined in (11171) . Finally the total Hamiltonian can be written 
as the sum of transverse part (11181) and longitudinal part (11211) 



, 2 



d^r^ ^ + 



D -P 



(V X A)^ 



2£o 2/i, 







where 



+ r du f d^k huBl{k,uj)Bxik,uj). 

A=l,2,3 

(122) 

^— r du I d^k y \F{uj)Bx(k,Loy^-'' + H.c]e),{k) (123) 

■n-Vz Jo J 



(2^)2 ^0 J ^^^ 2,3 

Therefore the minimal coupling Hamiltonian is obtained from the La- 
grangian ([1]) by using a canonical transformation then the two Hamiltonians 
in relation — (1301) and fl77|) are equivalent. 
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We can check the consistency of minimal couphng method with damped 
polarization method by comparting the obtained susceptibility from these 
two methods. In the damped polarization model the imaginary part of the 
susceptibility Imx{uj) is obtained as 



JmxH = ^^(aM + (124) 



where we have used Eq. ( 1561) and the definition of g{uj) in (l50l) . In minimal 
coupling method Imx{i^) can be written as 

Imxico) = (aiu) - /?M)^ (125) 

2puoeo 

where we have used fl89l) and flll7p . Now from fj48|) we can easily prove 

{a{u) - (3{u;)r = ^4(«H + /5(u;))^ (126) 



Substituting ffT26|) in ffT25|) . relation ffmj) will be obtained. 



4 Comparing the results of the minimal cou- 
pling method with other methods 

In this section we want to complete the equivalence of the minimal cou- 
pling method with other methods and show that the results obtained in the 
minimal coupling method, fl96|) — f fT02|) . are equivalent to those obtained in 
other methods. In equations fl96|) — f fT02|) the vector potential is obtained 
as a combination of the free field operators and the reservoir field operator 
respectively as 

A(r,t) = / + //, (127) 

where by setting the initial time at — oo, we can write the expressions for / 
and // as 



lim 




A=l,2 



2+(a;k, t - t')e*''"aA(k, t') + H.c.]ex{k), 



:i28) 
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and 



11= ,lim - / -7r\T / oo[aco,u^,t)B,ik,uj,ty'^-^ + H.c.]e,ik). 
t'^-oo So ^.^J (27r)2 Jo 

(129) 

In other words in Eqs. fl96l) and (110 II) we change time variable from t to 
t — t' and then let t' —oo. 

The explicit time dependence of z{u!k,t — t') and ^{(jj\^,uj^t — t') can be 
obtained from the the inverse Laplace transform formula 

fit) = L-'im] = hm l^^^^e^^-^^^'dco, (130) 
where f{t) is an arbitrary function. We have 



lim z^{u]f_, t — t') = 

^^o+t'^-ooy_oo 27r ujI + {iu + r]y[l + x{-^)] 

In equation (11311) the relation x(^^) = x(~'^) have been used where x(s) 
and x(co') are the Laplace and Fourier transforms of x(t) respectively. 

In equation (I13ip we have two situations: (i) The medium is nondispersive 
and X is a constant. In this case we use the limit t' — for simplicity and 
solve the equation (I13ip by using residue calculations and find 

z+{uj, t) = {l + n) e(^*) + (1 - n) e'^^*), (132) 

where n = ^/e which coincide with the result of [9]. 

(ii) The medium is dissipative. In this case according to the Kramers— Kronig 
relations x(a') is a complex function of u. Since t — t' > 0, the poles with 
positive imaginary part are important. Let there be N poles with Ui > 
then 



N 



lim z+{uk,t-t') = V lim an{uJk)e" 

n=l 
N 

= V lim [a„(a;fc)e^'^"'-(*-*')-'^"'*]e'^"^*', (133) 



, t'~*~oo 
n=l 
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where 

a^M = lim - -„){#^,FX^>- (134) 
From Eq. fll33p it is clear that 

hm z+{uj^,t-t') 0. (135) 

t'^—oo 

The same procedure can be used to obtain the exphcit form of the second 
term of equation (I127p and the final result is 



X lim lim 



,,^0+ t'^-oo_7^ 27r {i{uj + uj') + r]}{ujl - {u - iriY[l + xi.-^)]}' 

(136) 

In equation (I136p there is one real pole and the other poles have imaginary 
parts and tend to zero. Now by using the calculus of residues we obtain 



eoA^2^o y (27r)2 {ujl-uj'^e{uj)) 
X (5A(k,a;)e*(-*-'^'') + //.c.)eA(k), (137) 

where we have defined B\{\s.,uj) as 

lim SA(k,u;,t')e^'^*' = 5A(k,u;). (138) 

t'^—oo 

Using fl89|) to write F{uj) in terms of Imx{^) equation fll37p can be written 
as 



\87r^£o/ J J x^,2 k c2 _ tu"=£:(to') 
X {Bxik, a;)e^("*-'^ '-) + H.c.) exik). (139) 
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This equation is the same equation obtained in the phenomenological and 
the damped polarization model [51 E]. 

By the same method the longitudinal part of electric field can be obtained 

as 



87r%o/ Jo J 1 + xM 

X (53(k,a;)e^("'-)-/i.c.)e3(k), (140) 

„ : II 

which is in the phenomenological model and also the longitudinal part 
of the electric filed in the damped polarization model. 



5 Comparing the different methods 

In subsection 2.3 we showed that the phenomenological method can be ob- 
tained from minimal coupling method and in section 3 we obtained the 
minimal coupling Hamiltonian from the Lagrangian of damped polarization 
method. In addition as can be seen from (I55l)-(I59D, (^^-(^B> and (l96l)-(|T02ll. 
these three methods lead to the same results. Then in fact these three models 
are equivalent. In other word they are different techniques for solving the 
equation of motion of the same Lagrangian. 

Although these three methods have the same results and they are equiv- 
alent but each of them have its own advantage and disadvantage. 

The merit of damped polarization method with respect to the other meth- 
ods is that it is based on a Lagrangian and then we can use the standard 
canonical quantization method and the ETCR between different operators 
can be deduced from its Lagrangian. But by comparing the solution tech- 
nique of three methods we see that the solution technique of this model, 
which is based on Fano technique and diagonalization of the Hamiltonian, is 
hard and lengthy. In addition the extension of this model to nonisotropic and 
nonhomogeneous polarizable and magnetizable medium has not been down. 

The preference of phenomenological method is that its way of solving the 
equation is based on the Green function of the classical Maxwell equations. 
This correlation function would be useful in calculating the vacuum effect, 
like spontaneous emission, Casimir effect and Van der Wales force [T5] — p!8]. 
In addition the extension of this model to a nonhomogeneous and non- 
isotropic medium is easy and only we should solve the classical Green func- 
tion of the medium with nonhomogeneous and nonisotropic susceptibility [19] . 
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Also the extension of this model to a magnetizable medium is down by adding 
new noise operators for the loss of magnetization [2U]. But the disadvantage 
of this model is that it is not based on a Lagrangian and the commutation 
relation between vector potential and electric field is postulated and should 
be checked after quantization. In addition in [21] it is shown that in the min- 
imal coupling method the spontaneous emission can be solved easier than 
this model. 

The advantage of minimal coupling model is that in this model Maxwell 
equations are obtained from Heisenberg equation. Also, as mentioned in 
( !9T|) . the commutation relation between noise operators are obtained from 
the Heisenberg equation which give a better understanding of the nature of 
noise operators. 

Another merit of this method is that it can be extended easily to a nonho- 
mogeneous, non isotropic and nonlocal polarizable and magnetizable medium 
[9J and [22]. 

The disadvantage of this model is that since this model is not based 
on a Lagrangian then the commutation between vector potential A and the 
displacement vector D is postulated. 

Therefore up to the range of applicability of these three models, we could 
prove that these models are equivalent in a spatially homogeneous and non 
magnetized medium. 

6 conclusion 

In this paper we have shown that the minimal coupling method is equivalent 
to the Huttner-Barnet and phenomenological approaches up to a canonical 
transformation in a nonmagnetic medium. The magnetic properties of the 
medium are also included in the minimal coupling method contrary to the 
other methods. So for a general comparison, an extension of the Huttner- 
Barnet model to the case of a magnetodielectric medium is needed which is 
under consideration. 
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